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ABSTRACT 



Lumped-parameter and finite-element-models are devel- 
oped for the dynamic behavior of plane straight sided 
frames. The models do not include axial and shear defor- 
mation, but the equations of motion developed using the 
models allow for time varying external loading. The per- 
formance of these models is evaluated by a comparison with 
a standard transfer matrix method for the special case of 
free undamped vibration. The finite-element-model proves 
to be much the better model. For the first five modes, 
the finite-element-model with 27 degrees of freedom differs 
by no more than 0.4% from the "exact values" given by the 
transfer matrix method. The lumped-parameter-model with 
39 degrees of freedom gives errors roughly ten times as 
great. 
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NOMENCLAUTRE 



Cross sectional area of member i. 

Spectral matrix from diagonalization of J. 
Matrix relating d to c. 

Coefficient vector. 

Element of coefficient vector 

Modal matrix from diagonalization of J. 

Deflection and slope vector. 

Modulus of elasticity of member i. 

Flexural rigidity 
External force vector. 

Vector of x components of external force, 
x component of external force at node j. 
Vector of y components of external force, 
y component of external force at node j . 

d 2 w 

Row vector in powers of r for — . 



dr 



Row vector in powers of r for 4>, 

Identity matrix. 

Second moment of area A^ with respect to a centroidal 
axis perpendicular to plane of frame. 

Matrix sum of inertia terms. 

Inertia matrix. 

Mass moment of inertia of segment j about an axis 
normal to the plane of the structure and thru the 
segment centroid. 

Rotational inertia between j and j+1. 

Stiffness matrix. 
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kj = Stiffness of segment j. 

L = Incidence matrix. 

s . 

1 . = Length of segment of member i, = — . 

1 

M = Mass matrix. 

M = Rotational mass matrix. 

— r 

= Translational mass matrix, 
rru = Segment mass. 

= Moment vector. 

= Moment applied between node j and j-1. 

N = Total number of subdivisions. 

n. = Number of subdivisions in member i. 

1 

P = Matrix used to apply end conditions in lumped- 
parameter-model . 

= Row matrix in powers of r for w. 

Q = Moment arm matrix for external force vector. 

£ = Generalized force vector. 

R = Corner transformation matrix, 

r = Axial coordinate. 

_JL 

£ - Matrix product of B DB 2 . 

s^ = Length of member i. 

T = Kinetic energy. 

T^ = Rotational kinetic energy. 

T t = Translational kinetic energy. 

U = Potential energy. 

u = Vector of x components of displacement 
u = Axial displacement. 
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Uj = x displacement of node j . 
v = Vector of y components of displacement. 

Vj = y displacement of node j . 

= Specific weight of member i. 
w = Transverse deflection. 

X = Matrix of the x components of distance. 

X j = x-coordinate of the centroid of segment j . 

x^ = x component of distance from node j to node i+1. 
Y = Matrix of the y components of distance. 

Yj = y-coordinate of the centroid of segment j. 

= y component of distance from node j to node i+1. 

T t 

= product of B and D 






(x,y,z) = Cartesian coordinate system. 
a = Deviation angle at a corner of the frame. 

3 = cos a. 

Y = Angle between axial direction in local coordinate 
system and x direction in Cartesian system. 

n = sin a. 

0_ = Relative rotation vector. 

0^ = Relative rotation of hinge j. 

y = Mass per unit length. 

<p - Slope. 

a) = Natural circular frequency. 

Subscripts 

a = assembled matrices, 

i = Property of member i. 

j = Property at node j . 
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ij = Property from node i to node j in finite-element- 

model; from node j to node i+1 in lumped-parameter- 
model. 

u = Upper rows of a matrix as determined by end con- 
ditions . 

H = Remaining rows of matrix. 

Superscripts 

t = Transpose of vectors and matrices. 

" = Matrix after application of end conditions. 

* = Vector resulting from corner transformation. 

= In lumped-parameter-model denotes a once reduced 
matrix; in finite-element-model denotes the proper- 
ties of the element to the right of the corner. 
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INTRODUCTION 



When designing a structure, the engineer must be aware 
of the many dynamic disturbances that may act upon this 
structure. He must also be aware of the response that the 
structure will have to these disturbances. The calculation 
of the response is usually a very difficult undertaking. 

To assist in this calculation, the structure is generally 
characterized by a simplified model, and the analysis 
carried out using this model. 

In this paper, two models for analyzing plane frame 

structures with "chain" topology* are developed. The 

first model lumps both mass and flexibility and the second 

is a fi nite -element- mo del . To test these models, analyses 

to find the in-plane circular frequencies and translational 

displacements are carried out on a typical gable bent. The 

results of these analyses are compared to frequencies and 

displacements obtained through the use of the so-called 

1 ** 

transfer matrix method. 

The equations of motion as developed for the models 
provide for time varying external loads, but do not include 
damping. The analysis of the gable bent is carried out for 
the special case of free vibration. This simplification 



*Plane framed structures having 2 ends and a single 
path from end to end. 

**Superscripts refer to the bibliography which appears 
on page 48. 
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is required so that the principal modes obtained through 
the use of the models may be compared with values from the 
transfer matrix method. 
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DEVELOPMENT OF THE LUMPED-PARAMETER-MODEL 



In order to find in-plane frequencies and displace- 
ments of a structure consisting of straight sides with uni- 
form cross section, a model lumping both mass and flexibil- 
ity is developed. This model neglects axial and shear de- 
formation, but can include rotatory inertia. 

The model is formed by subdividing each straight mem- 
ber into segments. These segments are considered to have 
square ends. The length along the centerline of the 
straight member and the sum of the lengths of the segments 
replacing it are the same, even though this would cause some 
overlapping at the corners. The mass and flexibility of 
each segment is lumped at the segment centroid. The struc- 
ture is then replaced by a chain of straight rigid bars, 
pin connected at the locations of the segment centroids. 

The mass of the segment remains concentrated at the hinge, 
and the flexibility of the segment is replaced by a tor- 
sional spring. The equations of motion are developed and 
the model analyzed to obtain its natural frequencies and 
the displacements of the hinges. 

Rigid body translation is prevented by taking the left 
end as fixed or pinned. The right end is allowed to assume 
any of the common end conditions: fixed, pinned, simply 
supported, or free. 

Development of the Equations of Motion 

To illustrate the method, a structure containing two 
straight members intersecting at an acute angle is used. 
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Frame 

Figure 1 

Both ends of the structure are taken as fixed (Figure la) . 
For each straight member i (i=l,2), let 
= specific weight, lb/cu.in. 

= cross sectional area, sq.in. 

I. = second moment of the cross sectional area with 
respect to the z axis, in. 4 

= modulus of elasticity, lb/sq.in. 
s^ = length of member, in. 
n^ = number of segments; n^=2 , n 2=3. 
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Each section of the structure is divided into n^ seg- 
ments (Figure lb) . The n^ need not be the same for the dif- 
ferent members, but in each member the segments are of uniform 
length. The parameters of each segment are then lumped at its 
centroid. These parameters are 

= segment length = s^/n^ inches. 



2 

nu = segment mass = (AW£) ^/386 lb-sec. /in. 

= mass moment of inertia about an axis normal 
to the plane of the structure and through the 
segment centroid, 

= ir^uJ/12 + I i /A i ) in-lb-sec? 
k. = segment stiffness = (EI)./i. lb-in. 

1 Jr 1 

The structure is then replaced by the rigid bars and tor- 
sional springs (Figure lc) . The springs are unstrained in the 
datum position. The Cartesian coordinates of the end points 
and hinges of the structure in the datum position are taken 
as Xj and Y^ where j=l,2 ,3 ,4 ,5 ,6 , 7 . Thus the displacements 
of the hinges (u^, v^) in the positive x and y directions for 
a relative rotation, 0^, at hinge j are given by 
u = -Y 0 



( 1 ) 



v = X 0 



where 



u = 


col (U£ , 


u 3 ' 


u 4 ' 


U 5 ' 


V 


U 7 ) 


V = 


col (v 2 , 


v 3 ' 


v 4 ' 


v 5 ' 


v 6 ' 


V 


£ = 


col (0 1 , 


CD 

to 


0 3' 


CD 


V 


V 
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1 — 1 
>1 
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0 


0 


0 


0 




y 21 


y 22 
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0 
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0 




y 31 


y 32 


y 33 


0 


0 


0 


Y = 


y 41 


y 42 


y 43 


y 44 


0 


0 




y 51 


y 52 


y 53 


y 54 


y 55 


0 




[*61 


y 62 


y 63 


y 64 


y 65 


y 66_ 




X 11 


0 


0 


0 


0 


0 




X 21 


X 22 


0 


0 


0 


0 


X = 


X 31 


X 32 


X 33 


0 


0 


0 

t 

) 




X 41 


x 42 


X 43 


X 44 


0 


0 

i 




X 51 


X 52 


X 53 


X 54 


X 55 


0 




U|61 


X 62 


X 63 


X 64 


X 65 


X 6£ 


Y ij : 


= Y i + 1 " 


Y . 
D 








X . . : 

ID 


= X i + 1 - 


X . 
D 









The following notation has been used in this equation: 

1. Column vectors are denoted by lower case symbols. 

2. Rectangular matrices are denoted by capitalized symbols. 

3. All vectors and matrices are underlined. 

This notation will be used in all other such equations. 

Associated with each bar is a rotational inertia made up 
of half the mass moments of inertia of the preceding and fol- 
lowing segments. This rotational inertia is represented by 



18 



j . . = i(j . + j • . *i ) • 

J J 2 J 3+1 

Using equation (1), the kinetic energy, T, and the poten- 
tial energy, U, of the structure are formed. 



T = ie X MX6 + j-e Y MY6 + L J q L8 
U = 



(2) 



where 



M = diag(m 2 , m 3 , m 4 , , m g , 0) 

'L 0 = dia ^( J ii' J 22 ' J 33 ' d 44 ' J 55 ' J 66 ) 
K = diag (k^ , k 2 , k 3 , k 4 , k 5 , k g ) 



10 0 0 
110 0 
1110 
1111 
1111 
1111 



0 

0 

0 

0 

1 

1 



0 



0 

0 

0 

0 

1 



Derivatives with respect to time are denoted by dots and the 
transpose of a matrix by a superscript t. 

The Lagrangian equations of motion, 

d /3T\_ 9T 3U _ 

^ ( 3 ef *1 se - 2 

may be written 

J6_ + K6 = q (3) 

where 

J = X fc MX + Y^MY + 

q = generalized force vector. 
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Let the x components of the external forces applied at hinges 
2 to 7 be represented by the vector 



-- Col(f x2' f x3 ' f x4 ' f x5 ' f x 6 ' 



— x 



-x7' 



In a like manner define a vector _f . For the counterclockwise 
external moments applied between hinge i -1 and i, the vec- 

tor 



hi ~ col ( * 1^2 / K 3 / ^.4 > ^5 > ^g • 



is defined. The generalized force vector may be expressed in 
terms of these vectors as 



2 = L Y t f x + X fc f y . (4) 

It should be noted that the matrices J and K are both 
symmetric. This condition results when any symmetric matrix 
is postmultiplied by a matrix and is premultiplied by its 
transpose. This symmetry property is highly desirable when 
dealing with matrix eigenvalue problems and should be preserved 
if possible when applying boundary conditions. 

Application of Boundary Conditions 

Up to this point the development has been general in na- 
ture. We now restrict further development to free vibration. 
With this restriction, the only external forces applied are 
those required to satisfy end constraints. If the model has 
both ends fixed, the end conditions and corresponding exter- 
nal forces are 

Qj = 0, h\ 2 jt 0 

V 7 “ °> f y 7 * 0 

u 7 = °' f x7 * 0 

®1 + ®2 + ®3 + ^4 + + ®g = O' ^7 ^ 0 • 
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Equation (4) can therefore be reduced to 



where 



f = 



Q = 



g. = Q f 



col (h\ 2 , 


f -t t 

y 7 


fl 


X 61 


*61 


i 

0 


X 62 


^62 


0 


X 6 3 


^63 


0 


X 64 


^64 


0 


X 65 


^65 


0 


X 66 


^66 



f x7' V 

1 

1 

1 

1 

i 

1 : 

* 

i 

l 



( 5 ) 



The columns of Q are the first and last columns of L*" , and 
the last columns of and Y t . 

The end conditions may be represented by 




( 6 ) 



This equation is derived from the following set of equations. 



e 1 = 0 



u 7 “ X 6i e i +X 62 0 2 +X 63 0 3 +X 64 e 4 +X 65 0 5 +X 66 0 6 = 
v 7 = y 6lV y 62 e 2 +y 63 e 3 +y 64 e 4 +y 65V y 66 0 6 = 



0^ + 02 + 63 + 0^ + 0^ + 0g = 0 



0 

0 



A process must be devised to incorporate the boundary 
conditions. This process will reduce the number of degrees of 
freedom by an amount equal to the number of constraint condi- 
tions, in this case four. The result will be a pair of homo- 
geneous equations with two remaining 0j's. 
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A normal Gaussian elimination technique for the solution 
of simultaneous equations is used to solve for f^ , f ^ , I , 
and . Through the use of the end conditions, the number of 
unknown 0^ is next reduced to two. The four non-homogeneous 
equations containing the forces and moments can be deleted, 
leaving only the two homogeneous equations. A similar elim- 
ination could be used for any other set of end conditions. 

This process can be represented by a set of matrix opera- 
tions. Equation (3) is written in partitioned form as 



J 




r 

K 


i 

j 


— i 

Q 


— u 


0 + 


— u 


£ = 


— u 

i 

i 


1 

■dI 

1 




A 




2* 



where the subscript u represents the rows which will contain 
the non-homogeneous equations, and £ the remaining rows. This 
may require a row interchange to insure a non-singular Q . 

In the present case is 4x4. Gaussian elimination is now 
performed to solve for the forces and moments. In this proc- 
ess the elements of will be reduced to zero. The lower 
portion of equation (7) may be rewritten 



iii* + = 0 (8) 

where the prime denotes reduced matrices. 

The end conditions, equation (6) , may be rewritten 



and from this we get 

Quia + - 0 



t t 

Q Q „ 

— u — £ 
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and therefore 



6 = - (Q t ) V 0 o 

— u — u —i—i 



Thus 0_ u may be eliminated from _0 by the transformation 



e = 






h m 



(9) 



where is the identity matrix, in this case of order 2x2. 
It is easily shown that 



P fc Q = 0 

and that in equation (8) 



and 

■si = ■ 

The 0_ may then be replaced in equation (8) by P6 „ . After dif- 
ferentiating equation (9) twice, the £ of equation (8) may 
also be replaced, and the following expression obtained: 

+ - 0 (10) 

where 

Jj = E^JP 
= P fc Kp . 

The boundary conditions can therefore be applied and a set of 
homogeneous equations obtained by this process. This trans- 
formation does not destroy the symmetry of the matrices. 
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Eigenvalue Solution of the Equations of Motion 

Assuming that the vibration will be sinusoidal, equation 
(10) may be written 

< n > 

where co is the natural circular frequency. This equation is 

in a form suitable for eigenvalue solution. 

The solution of equation (11) is now undertaken. This 

2 

solution will involve eigenvalues and eigenvectors. A sym- 
metric matrix may be expressed in terms of its modal matrix 
and its spectral matrix where the spectral matrix is a diag- 
onal matrix of eigenvalues and the columns of the modal matrix 
are eigenvectors. This property is used several times in the 
solution. A further important property of the modal matrix 
is that the eigenvectors are normalized with respect to the 
identity matrix. Thus the transpose of the modal matrix is 
also its inverse. 

If is resolved into its (diagonal) spectral matrix 
B and modal matrix D according to the relation 

Jj = DBD t , 

equation (11) becomes 

= 0J 2 DBD t e jl . 

Since DD*" = I, the above equation may be rewritten 

KpD t 9 „ = gj 2 DBD t e „ . 

Since B is a diagonal matrix of positive elements , we define 
the principal square foot as B2 where the elements of this ma- 
trix are the positive square roots of the elements of B. Not- 
_i_ i 

ing that B 2 B 2 = 1 , the following form results, 
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-lit 2 1 1 t 

K£DB 2 B2D l '_e £ = u DB2B2D 0^ . 

, 1 t 

Letting z = B^D 0^ , 

_1 9 I 

K^DB- 2z = co DB2z . 

2 

Solving for co z, 

l +- -1 2 

B 2D KJDB 2z = co z . 

_1 t _i 

Letting S_ = B 2D K^DB 2 f the equation may be rewritten 

2 

Sz = co z_ . 

2 

By diagonalization of S, the values of co and z_ may be ob- 
tained. From these values, the natural frequencies are ob- 
tained directly. The values of the rotations of the nodes 

_i 

are obtained from the equation 6^ = DB2£. The rotations 
at all the hinges of the structure can be recovered by use 
of equation (9) . The translation displacements of the hinges 
are easily obtained with the use of equations (1). 
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DEVELOPMENT OF FINITE-ELEMENT-MODEL 



General 

A finite-element- model is developed to describe the 
in-plane motion of a frame with straight members of uniform 
cross section. The model neglects axial and shear deforma- 
tion and can include rotatory inertia if desired. 

The use of a finite element approach requires the divi- 
sion of the structure into elements* and the determination 
of element mass and stiffness matrices. The determination 
of the elemental mass and stiffness matrices will be taken 
up first. 

Consider a general element of length i, in a local 
coordinate system. In this system, r is the axial coor- 
dinate, w is the transverse deflection, and <f> is the slope. 
The element has no flexibility in the r direction. If we 
assume that w may be approximated by a third degree polyno- 
mial in r, then 

2 3 

w = c + c.r + c~r + c-,r 
O 1 Z J 

and 

<j> = = c. + 2c„r + 3c.>x 2 

r dr 1 2 3 

where c q , c^, c 2 , c 3 are unknown coefficients. These equa- 
tions may be expressed as 

t 

W = £ C 

cj) = h^c 

where 

= row(l, r, r 2 , r 3 ) 
h fc = row(0, 1, 2r, 3r 2 ) 
c = col(c Q , c^, c 2 , c^) . 
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If w and <p are evaluated at the ends of the element, the re- 
sults may be expressed as 

d = 

where 



d = col(w 0 , <j> 0 , w £ , <}> £ ) 



-1 



1 

0 

1 

0 



0 

1 

£ 



0 0 
0 0 



a 2 £ 3 



1 2 £ 3 £ 2 



and the subscripts denote evaluation at r= 0 , and r=£. By 
inverting C ^ we may solve for c, 



c = Cd 



( 12 ) 



where 



£ 1zl Zi 



; £ £ 

2 1 



0 0 
3 -1 



-2 



£ 3 £ 2 £ 3 



1 

o 2 



Bending Energy 

5 The bending energy, U, of the element is given by 



, A A2 \ 2 

u = 

2 'o \dr“ 



dr 



( 13 ) 



where El is the flexural rigidity for in-plane bending. 
Now 



d 2 w 0 , , 

— 7 ? = 2c~ + 6c 0 r = t 
2 2 3 g c 



dr 



where 



g fc = row(0, 0, 2, 6r) . 
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